arXiv:1509.00787v2 [physics.space-ph] 18 Dec 2015 


Quantifying Gyrotropy in Magnetic Reconnection 

M. Swisdalfl 

A new scalar measure of the gyrotropy of a pressure tensor is defined. Previously suggested 
measures are shown to be incomplete by means of examples for which they give unphysical results. 
To demonstrate its usefulness as an indicator of magnetic topology, the new measure is calculated 
for electron data taken from numerical simulations of magnetic reconnection, shown to peak at 
separatrices and X-points, and compared to the other measures. The new diagnostic has potential 
uses in analyzing spacecraft observations and so a method for calculating it from measurements 
performed in an arbitrary coordinate system is derived. 


I. INTRODUCTION 

The term magnetic reconnection, as commonly under¬ 
stood, refers to a process whereby the energy of the em¬ 
bedded field is transferred to the surrounding plasma 
via bulk acceleration and Ohmic dissipation. While fre¬ 
quently associated with changes in the topology of the 
magnetic field, precisely defining where and how recon¬ 
nection occurs is a surprisingly subtle task [T]. 

In what are known as 2.5D geometries, where varia¬ 
tions in one direction are suppressed, the coupling of the 
reduced dimensionality with the divergence-free nature of 
the magnetic field allows field lines to be completely char¬ 
acterized as contours of a scalar function, ip. {ip is equiv¬ 
alent to the component of the magnetic vector potential 
in the invariant direction.) Saddle points of ip, known 
as X-points, mark locations of both topological change 
and magnetic reconnection. In fully three-dimensional 
systems, field lines can no longer be described by a single 
function and characterizing reconnection becomes signif¬ 
icantly more difficult. Finding the locations where mag¬ 
netic reconnection occurs is non-trivial, even in numerical 
simulations mg. 

Yet simulations have the luxury of a synoptic view 
with, in principle, access to the complete description of 
the plasma at any point in the domain. Spacecraft ob¬ 
servations, in contrast, must make due with a limited 
set of measurements available from, at most, a few loca¬ 
tions. The topological identification of reconnection sites 
is non-local in nature and hence extremely difficult to 
implement with these restrictions. Determinations must 
instead be based on local measurements of the electro¬ 
magnetic fields and particles. 

Due to their small Larmor radii and brief gyroperiods, 
electrons are closely tied to the magnetic field, much more 
so than the heavier ions. By streaming for long distances 
along field lines they efficiently probe magnetic structure 
and so locations of topological interest, such as recon¬ 
nection sites, should leave signatures in electron mea¬ 
surements. [g identified one such signature by noting 
that, because of finite Larmor radius effects, the guiding 
center approximation breaks down near an X-point. One 
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consequence is that reconnection, at least in 2.5D, re¬ 
quires an asymmetry (specifically, non-gyrotropy) in the 
electron pressure tensor there. Spacecraft measurements 
of electron distribution functions can hence be used as 
a proxy for topology. In this paper we present a diag¬ 
nostic quantifying the size of departures from gyrotropy 
and demonstrate, through numerical simulations, that it 
faithfully maps regions associated with magnetic recon¬ 
nection. 


II. MEASURING GYROTROPY 


Given a distribution function /(x, v) that describes the 
density of particles with mass m in a phase space defined 
by position (x) and velocity (v), the pressure P can be 
defined as 

P = m / ww/d^v , (1) 


where w = v — v is the velocity relative to the mean flow 
V. While P is, in general, a tensor with six independent 
components, it can be characterized by a single scalar 
when / is spherically symmetric (e.g., a Maxwellian). A 
magnetic field, by imposing a preferred direction, can 
produce a cylindrically symmetric, or gyrotropic, distri¬ 
bution. In this case the pressure tensor takes the form 

P = Pj|bb-kPj.(I-bb) , (2) 

where b = |B|/P is a unit vector in the direction of the 
magnetic field, Py and Pj_ are scalars, and I is the unit 
tensor. In the equivalent matrix formulation. 


P = 



0 0 \ 
p± 0 
0 P± J 


( 3 ) 


Observationally, pressure tensors are never completely 
gyrotropic because of the presence of non-zero off- 
diagonal components, although these are often much 
smaller than the diagonal terms. Since departures from 
gyrotropy are expected to be associated with locations 
of magnetic reconnection, a quantitative measure of this 
smallness would be useful. 

With this goal in mind, [6] proposed a measure, dubbed 
the agyrotropy and denoted A0e, that characterizes a 
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distribution function’s weighted average of dispersions 
of velocities perpendicular to the local field direction 
(see Appendix A of their paper for a full description). 
It varies between 0 and 2 with the minimum suppos¬ 
edly corresponding to gyrotropic distributions. However, 
as demonstrated below, there exist agyrotropic pressure 
tensors for which A 0 f. = 0. 

[7] proposed another measure, denoted D^g and named 
non-gyrotropy. It is proportional to the root-mean- 
square of the off-diagonal elements of the pressure tensor 
normalized by the local thermal energy. As with the agy- 
rotropy, D^g = 0 for a gyrotropic distribution but, as also 
shown below, there exist pressure tensors with maximal 
departures from gyrotropy for which Dng is arbitrarily 
close to 0. 

A brief mathematical digression is necessary in order to 
justify a better measure. Because pressure tensors orig¬ 
inate as moments of distribution functions, their matrix 
representations must satisfy certain constraints. Specif¬ 
ically, a physically meaningful pressure tensor is a real 
3x3 symmetric matrix (i.e., Pij = Pji) with non-negative 
eigenvalues. Matrices possessing these properties, includ¬ 
ing all pressure tensors, are called symmetric positive 
semi-definite. If all of the eigenvalues are positive the 
matrix is positive dehnite. (That the eigenvalues of a 
pressure tensor must be non-negative follows from the 
requirement that after diagonalization, which is always 
possible for a real symmetric matrix, the matrix elements 
are both the eigenvalues and the components of the pres¬ 
sure along the basis vectors of the rotated frame. For 
physical distributions the latter, and hence the former, 
must be greater than or equal to zero.) 

Defining a scalar measure of gyrotropy requires a re¬ 
sult from the theory of symmetric positive semi-definite 
matrices. It can be shown [8] that a matrix is positive 
semi-definite if and only if all principal minors - the de¬ 
terminants of every submatrix arising from deleting the 
same set of rows and columns - are non-negative. (If zero 
eigenvalues are excluded, so that the matrix is positive 
definite, Sylvester’s criterion applies [S]. In this case, 
a necessary and sufficient condition is that the leading 
principal minors, i.e., the determinants of all upper-left 
submatrices, be positive.) 

A pressure tensor with gyrotropic diagonal, but non¬ 
zero off-diagonal, entries can be written as 


= P 


4^11 

P12 

Pl 3 

P12 

P± 

P23 

Pl 3 

P23 

P± 


(4) 


In this representation one of the coordinate axes points 
in the direction of the local magnetic field and the others 
are oriented such that the final two components of the 
diagonal of P are equal. (That this is always possible is 
demonstrated in the Appendix.) Since P is positive semi- 
definite the result discussed above implies the inequalities 

P?2<P\\P± P?3<P\\P± Pis < Pi (5) 


The converse does not hold: satisfaction of these inequal¬ 
ities does not guarantee that P is positive semi-definite. 

Adding the inequalities implies a natural definition for 
a measure of gyrotropy: 


Q 


P?2+P?3 + Pi3 

Pi + 2P_lP|| 


( 6 ) 


For gyrotropic tensors (5 = 0, while for maximal depar¬ 
tures from gyrotropy Q = 1. Evaluating Q from equation 
1^ requires P to be in the form given by equation even 
though, in general, P is measured in an arbitrary coor¬ 
dinate system. While it is always possible to rotate P 
into a frame in which it has the form of equation the 
Appendix demonstrates how to calculate Q without such 
a transformation by using tensor invariants. 

The theory of positive semi-definite matrices provides 
a firm mathematical basis for the dehnition of Q. With 
it in mind, it is possible to construct pressure tensors for 
which other proposed measures give physically unreason¬ 
able results. Consider a system with a pressure tensor, 
of the type shown in equation given by 

/I 1/2 0\ 

P = 1/2 1 0 . (7) 

\ 0 0 ij 


Since P has eigenvalues of 1/2, 1, and 3/2, it is positive 
semi-definite and hence physically valid. Yet, despite the 
off-diagonal elements that make it clearly not gyrotropic, 

A 0 , = O. (8) 


For reference, ^/Q = a/S/G and Dng = V2/3 in this case. 
(The other gyrotropy measures considered here scale as 
ratios of first powers of pressure components while Q 
scales quadratically. The subsequent discussion, where 
necessary to ensure an unbiased comparison, uses y/Q.) 

In fact A0e = 0 for any positive semi-definite P in 
the form of equation with P 23 = 0, irrespective of the 
other off-diagonal components. This occurs because, by 
construction, A 0 e only measures departures from sym¬ 
metry in the plane perpendicular to the magnetic field 
axis. However, as equation demonstrates, there are 
other ways in which pressure tensors can depart from 
gyrotropy. The vanishing of A 0 e for such cases indi¬ 
cates that it does not fully describe a pressure tensor’s 
gyrotropy. 

The nongyrotropy measure proposed by [7] is defined 
as 


Dng — 


VlP?2+P?3+Pi3) 

P\\ + 2 P± 


(9) 


The most significant difference between it and ^/Q occurs 
in the normalization, which for Dng is to the thermal 
energy, P± +P^^/2. As with the other measures, Dng = 0 
for gyrotropic tensors. The maximum value of Dng varies 
with P\\/P±, which is not the case for either Q or A 0 . 

Its global maximum of occurs when Py = P_l. 
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The implications of this variation become apparent 
through consideration of the pressure tensor 

( X y/x y/x\ 

1 1 ( 10 ) 

1 1 / 

with a; > 0. Eigenvalues of 0, 0, and x + 2 mean that P 
is positive semi-definite. The gyrotropy measures are 

Dng = ^ ^ ^ ■ 

( 11 ) 

This pressure tensor is maximally non-gyrotropic for any 
value of X because any perturbation that increases the 
off-diagonal components will drive one of the eigenvalues 
negative. Any measure of the gyrotropy for this P should 
thus be independent of x and should take its maximal 
value. Yet Dng varies with x. This behavior does not 
depend on P having zero eigenvalues, although such a 
case is the simplest to analyze. Similar pressure tensors 
exist that both have three positive eigenvalues exist and 
approach arbitrarily close to maximal non-gyrotropy yet 
for which Dng varies in a manner similar to Equation |11[ 
The examples of equations and are instances 
where the dehciencies of A0e and Dng in quantifying 
gyrotropy are particularly clear. To a lesser or greater 
degree similar issues will be present for other pressure 
tensors, although in any specific example it may be dif¬ 
ficult to tease out the degree to which the gyrotropy is 
mischaracterized. 


III. SIMULATIONS 

While the arguments of the previous section suggest 
that Q is a good measure of gyrotropy, it is a separate 
question as to whether it is also useful as a proxy for the 
identihcation of reconnection sites. As a first attempt 
at an answer, we investigate particle-in-cell simulations 
performed with the code p3d [9] . It employs units based 
on a held strength Bq and density uq, with lengths nor¬ 
malized to the ion inertial length di = cjujpi, where u>pi is 
the ion plasma frequency, and times to the ion cyclotron 
time 


A. Run 1: Anti-Parallel Reconnection 

We hrst consider a variation of the standard GEM 
Challenge m- The simulation domain has dimensions 
{Lx,Ly) = (51.2,25.6) with an initial magnetic held of 
Bx = tanh(y/r(;o) and wq = 0.5. There is no initial 
held component in the z direction. To ensure pres¬ 
sure balance, the density n = ni, + sech^{y/ wq) , with 
rib = 0.2, and the electron and ion temperatures are ini¬ 
tially isotropic with = 1/12 and Ti = 5/12. The 
ion-to-electron mass ratio is 400 and the speed of light is 
40. The spatial grid has resolution Aa; = 1/160, which 


(A) 
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FIG. 1. Images of \/Q (panel A) and A0e/2 (panel B) 
with superimposed magnetic field lines (dotted) for the run 
described in section |III A| Colors represent different values 
in each image; the central bar shows the relative variation. 
The data have been averaged over 0.03f2~^ = 1212^^, or ~ 2 
electron Larmor orbits in the asymptotic field, and over 0.5de 
in each direction. The panels on the right show cuts along 
the corresponding colored lines in the main images. 


means there are 8 gridpoints per electron inertial length 
{de) and « 3 per electron Larmor radius in the maximum 
field. There are 100 particles per cell for each species in 
the asymptotic region. 


Figure shows a comparison of \/Q and A0el‘2- dur¬ 
ing a time of steady reconnection {Dng closely follows 
a/Q in this case and so is not plotted). The magnetic 
field and pressure tensor components were averaged over 
« 2 electron Larmor orbits and 4 gridpoints (O.Sde) in 
each direc tion before computing the gyrotropy measures. 
Equation A8 was used to calculate y/Q-, A0e was cal¬ 
culated from the formulas given in Appendix A of [S]. 
In order to minimize subjectivity, no data manipulation 
(e.g., noise rejection), other than the spatial and tem¬ 
poral averaging, was performed. Only a fraction of the 
simulation domain is pictured in order to highlight the 
regions exhibiting strong departures from gyrotropy. On 
the right-hand side of each image are the traces of two 
cuts, one through the X-point (blue) and the other « 5di 
downstream (red). 

Both A0e and Q are significant near the X-point, al¬ 
though there are minor differences in the aspect ratio and 
structural details of the bright regions. The bifurcation 
at y/di — 0 and —2 < x/di < 2 mirrors that seen in the 
charge and current densities (not shown) and reflects the 
complicated dynamics of Speiser particle orbits near the 
region of null field. The values within a few gridpoints of 
{x, y) = (0,0) carry somewhat more uncertainty because 
gyrotropy is an ill-defined quantity when B = 0. 


While ^/Q and A0e/2 both vary between 0 and 1, the 
more meaningful comparison is between the asymptotic 
value of each cut, which gives a reasonable approximation 
of the inherent particle noise, and the peak value at the 
X-point. These are roughly equal for the two measures. 
On the other hand, the magnetic separatrices (red peaks) 
are comparatively much weaker in A0e and stronger in 
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FIG. 2. Images of y/Q (panel A) and A0e/2 (panel B), in 
the sa me format as Figure for the run described in section 


IIIB 


The data have been averaged over 0.511^./ = 12.5flc 


stream magnetic island while A0e quickly peters out. 
However, the strongest difference is in the background 
noise level. While the two measures have similar peak- 
to-background ratios for the simulation of Section |III A[ 
here they differ by a factor of « 4, with ^/Q smoother 
than in Figure [l] and A0f, noisier. The improvement in 
^/Q can be attributed to the strong guide field and its 
driving of the plasma towards gyrotropy. The indenta¬ 
tion in the magnetic field lines at x/di « 11 is a transient 
feature due to the interaction of the outflowing plasma 
and the remains of a small magnetic island expelled from 
the X-point earlier in the simulation. 


IV. DISCUSSION 


y/Q. The separatrices mark topological boundaries with, 
in principle, zero width. If the associated particle signa¬ 
tures exhibit scales near the Larmor radius (as is known 
to be the case around the X-line) then the intermingling 
between the upstream and downstream plasmas should 
produce departures from gyrotropy there. 

Both measures show humps just upstream of the X- 
point that indicate small departures from gyrotropy. 
These may arise from decreases in P± (due to mag¬ 
netic moment conservation) followed by scattering to 
other components. Another possibility is that beams 
of inflowing electrons, when superimposed on the base 
Maxwellian, drive the total distribution function away 
from cylindrical symmetry. 


B. Run 2: Guide-Field Reconnection 

As a second case we consider a force-free equilibrium 
where the initial density (n = 1) and temperatures 
{Ti = Tf, = 1/8) have no spatial variations. The re¬ 
connecting field has the form = tanh(?//r(;o), here 
with Wo = 1, but now the out-of-plane (guide) compo¬ 
nent Bz varies so that B^ + B^ is constant. The ini¬ 
tial guide field is asymptotically 2 and rises to y/E at 
the current sheet’s center. The domain has dimensions 
{Lx,Ly) = (51.2,25.6), the ion-to-electron mass ratio is 
25, and the speed of light is 15. The spatial grid has reso¬ 
lution Acc = 1/40, which implies 8 gridpoints per de and 
« 2 per electron Larmor radius in the maximum field. 
There are 100 particles per cell for each species. 

Figure in the same format as Figure shows the 
region near the X-point. Both ^/Q and A0f. peak there, 
although again the morphologies in the downstream re¬ 
gion exhibit differences. The asymmetry between the 
separatrices, a common feature of guide-field reconnec¬ 
tion, is apparent. As in Figurethe cuts show that the 
separatrices are somewhat stronger in than in A0f.. 
In images of the entire domain (not shown), the enhance¬ 
ment of ^/Q on the upper-right and lower-left separatrices 
continues around nearly the entire exterior of the down- 


We have proposed a new, mathematically rigorous 
measure of the gyrotropy of an arbitrary pressure tensor 
and presented numerical simulations of 2.5D reconnec¬ 
tion demonstrating that gyrotropy violations peak near 
magnetic topological boundaries. Testing the efficacy of 
Q in three-dimensional simulations and with spacecraft 
data are obvious next steps. It is unlikely that any sin¬ 
gle diagnostic can unfailingly and unambiguously iden¬ 
tify reconnection sites in spacecraft data, but Q offers 
a complementary approach to methods based on other 
measurements. The calculation of Q requires the full 
electron pressure tensor - available, for instance, on the 
recently launched Magnetospheric Multiscale Mission - 
and can be easily performed in any coordinate system 
using the algorithm presented in Appendix fA} 

While Q has some similarities to previously proposed 
measures, particularly the non-gyrotropy parameter D^g 
of [7], there are pressure tensors (e.g., equation 10) for 
which Dng does not correctly characterize the gyrotropy. 
However, in many instances the two give similar results. 
Their ratio can be written as 


Q _ {x-\- 2)^ 

^ " 8(1 + 2 x) 


( 12 ) 


where x = P||/Fj_. The right-hand side takes its min¬ 
imal value for x = 1 and only varies by « 30% for 
0 < a; < 4. Although the ratio can become arbitrar¬ 
ily large when a: —> oo, approaching that limit in a real 
system will trigger other effects. Of particular note is 
the firehose mode, for which the linear instability crite¬ 
rion is /3_l(P||/Pj_ — 1) > 2. Under normal circumstances 
plasmas cannot significantly exceed this threshold, which 
thus implies x < 1 + 2 /As a consequence, the dif¬ 
ferences between Dng and a/Q are starkest for /3_l <C 1. 
(This argument is strictly only applicable for gyrotropic 
plasmas because they are assumed in deriving the fire¬ 
hose stability boundary. It should be a reasonable ap¬ 
proximation when the pressure tensor is near gyrotropy, 
but significant off-diagonal components will likely modify 
the instability threshold.) 

The differences between Q and the agyrotropy A0f, 
are more significant. The latter only quantifies one type 
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of departure from gyrotropy, the breaking of cylindrical 
symmetry in the plane perpendicular to the magnetic 
field, while the former handles the general case. Of more 
pressing interest to the interpretation of simulations and 
observations, however, is how well the measures iden¬ 
tify regions of interesting magnetic topology. Both suc¬ 
cessfully illuminate X-lines for the simulations presented 
here, but Q does a better job of tracing magnetic sepa- 
ratrices, particularly in guide-field reconnection. Further 
analysis is necessary to make these comparisons quanti¬ 
tative. 

Finally, an interesting connection exists between P and 
the inertia tensor of a rigid body. Since both are sym¬ 
metric and positive semi-definite, a measure analogous 
to Q also exists for the former. To be relevant, however, 
something external to the body must define a preferred 
direction in a manner similar to the role the magnetic 
field plays in defining gyrotropy. 


Calculating K^PM and setting the second and third di¬ 
agonal components equal gives the required angle: 

tan(20) = (A4) 

In the new frame the final two diagonal components are 
{P±i + Fj. 2)/2. The off-diagonal elements must also be 
known in the new coordinate system in order to com¬ 
pute Q. Although explicit expressions for them can be 
calculated, the algebra is somewhat tedious. 

An easier method relies on the fact that certain combi¬ 
nations of tensor elements are invariant under coordinate 
rotations. Through their use, Q can be computed when P 
is known in any coordinate system. In particular, a 3 x 3 
tensor with elements Pij in an arbitrary Cartesian sys¬ 
tem with coordinates (x, y, z) has three invariants, two 
of which are the trace 
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Appendix A: Calculating Q 


and the sum of principal minors 


^2 — Pxx Pyy 


- PxxPzz + PyyPzz 
i^PxyPyx P PxzPzx 


PyzPzy). (A6) 


Beginning with a pressure tensor in an arbitrary co¬ 
ordinate system, the calculation of Q from equation 
requires a transformation into a frame in which the di¬ 
agonal components are in gyrotropic form. A tensor A 
transforms under coordinate rotations according to the 
prescription 


For symmetric matrices Pij = Pji and the final three 
terms of I 2 simplify. The third invariant, not used here, 
is the determinant. 

With these definitions Q can be expressed in terms of 
Ii, I 2 , and the parallel pressure Py. The latter can be 
calculated from 


A' = M^AM 


(Al) 


where K is a rotation matrix. Both the symmetry of 
A and its eigenvalues are preserved after rotation of the 
coordinate axes. 

As an example, we demonstrate that, beginning with 
a pressure tensor of the form 


^11 

Pa 

Pb \ 


Pa 

Pj.1 

Pc , 

(A2) 

Pb 

Pc 

Pj.2 / 



it is always possible to find a coordinate system in which 
the final two diagonal entries are equal. To do so, rotate 
around the magnetic field direction with the rotation ma¬ 
trix 


/I 0 0 \ 

M = 0 cos 9 — sin 0 . (A3) 

y 0 sin 0 cos 9 j 


P|l = b • P • b 

= b^Pxx + i'yPyy + b^PzzP (A7) 

‘^{pxbyPxy F bxh^Pxz P byh^Py^^ 


where bi is the component of the unit vector aligned 
with the local magnetic field. After some algebraic ma¬ 
nipulations. 


0 = 1 - 


(^1 


4/2 

^11)^1+3Pi|) ■ 


(AS) 


With this formula Q can be simply computed in an arbi¬ 
trarily oriented coordinate system, without need for co¬ 
ordinate transformations and the associated matrix mul¬ 
tiplications. 
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